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Abstract — Recently,  cooperative  communications,  in  the  form 
of  having  each  node  equipped  with  a  single  antenna  and  exploit 
spatial  diversity  via  some  relay  node’s  antenna,  is  shown  to 
be  a  promising  approach  to  increase  data  rates  in  wireless 
networks.  Under  this  communication  paradigm,  the  choice  of 
a  relay  node  (among  a  set  of  available  relay  nodes)  is  critical 
in  the  overall  network  performance.  In  this  paper,  we  study  the 
relay  node  assignment  problem  in  a  cooperative  ad  hoc  network 
environment,  where  multiple  source-destination  pairs  compete 
for  the  same  pool  of  relay  nodes  in  the  network.  Our  objective  is 
to  assign  the  available  relay  nodes  to  different  source-destination 
pairs  so  as  to  maximize  the  minimum  data  rate  among  all  pairs. 
The  main  contribution  of  this  paper  is  the  development  of  an 
optimal  polynomial  time  algorithm,  called  ORA,  that  achieves 
this  objective.  A  novel  idea  in  this  algorithm  is  a  “linear  marking” 
mechanism,  which  maintains  linear  complexity  of  each  iteration. 
We  give  a  formal  proof  of  optimality  for  ORA  and  use  numerical 
results  to  demonstrate  its  capability. 

Index  Terms — Cooperative  communications,  relay  node  assign¬ 
ment,  achievable  rate,  ad  hoc  network,  optimization. 

I.  Introduction 

SPATIAL  diversity,  in  the  form  of  employing  multiple 
transceiver  antennas,  is  shown  to  be  very  effective  in 
coping  fading  in  wireless  channel.  However,  equipping  a  wire¬ 
less  node  with  multiple  antennas  may  not  be  practical,  as  the 
footprint  of  multiple  antennas  may  not  fit  on  a  wireless  node 
(particularly  on  a  handheld  wireless  device).  To  achieve  spatial 
diversity  without  requiring  multiple  transceiver  antennas  on 
the  same  node,  the  so-called  cooperative  communications  has 
been  introduced  [10],  [16],  [17].  Under  cooperative  communi¬ 
cations,  each  node  is  equipped  with  only  a  single  transceiver 
and  spatial  diversity  is  achieved  by  exploiting  the  antenna  on 
another  (cooperative)  node  in  the  network. 

We  consider  two  categories  of  cooperative  communications, 
namely,  amplify -and-forward  (AF)  and  decode-and-forward 
(DF)  [10].  Under  AF,  the  cooperative  relay  node  amplifies  the 
signal  received  from  the  information  source  before  forwarding 
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it  to  the  destination  node.  Under  DF,  the  cooperative  relay 
node  decodes  the  received  signal,  and  re-encodes  it  before 
forwarding  it  to  the  destination  node.  Regardless  of  AF  or 
DF,  the  choice  of  a  relay  node  plays  a  critical  role  in  the 
performance  of  cooperative  communications  [1],  [2],  [24].  As 
we  shall  see  in  Section  III,  an  improperly  chosen  relay  node 
may  offer  a  smaller  data  rate  for  a  source-destination  pair  than 
that  under  direct  transmission. 

In  this  paper,  we  study  the  relay  node  assignment  problem 
in  a  cooperative  ad  hoc  network  environment.  Specifically, 
we  consider  an  ad  hoc  network  where  there  are  multiple 
active  source-destination  pairs  and  the  remaining  nodes  can  be 
exploited  as  relay  nodes.  We  want  to  determine  the  optimal 
assignment  of  relay  nodes  to  the  source-destination  pairs  so  as 
to  maximize  the  minimum  data  rate  among  all  pairs.  Although 
solution  to  this  problem  can  be  found  via  exhaustive  search 
(among  all  possible  relay  node  assignments),  the  complexity 
is  exponential.  Our  goal  in  this  paper  is  to  find  an  algorithm 
with  polynomial-time  complexity  to  solve  this  problem. 

A.  Main  Contributions 

In  this  paper,  we  study  how  to  assign  a  set  of  relay  nodes 
to  a  set  of  source-destination  pairs  so  as  to  maximize  the 
minimum  achievable  data  rate  among  all  the  pairs.  The  main 
contributions  of  this  paper  are  the  following. 

«  We  develop  an  algorithm,  called  Optimal  Relay  Assign¬ 
ment  (ORA)  algorithm,  to  solve  the  relay  node  assign¬ 
ment  problem.  A  novel  idea  in  ORA  is  a  “linear  marking” 
mechanism,  which  is  able  to  offer  a  linear  complexity  at 
each  iteration.  Due  to  this  mechanism,  ORA  is  able  to 
achieve  polynomial  time  complexity. 

•  We  offer  a  formal  proof  of  optimality  for  the  ORA 
algorithm.  The  proof  is  based  on  contradiction  and  hinges 
on  a  clever  recursive  trace-back  of  source  nodes  and  relay 
nodes  in  the  solution  by  ORA  and  another  hypothesized 
better  solution. 

«  We  show  a  number  of  nice  properties  associated  with 
ORA.  These  include:  (i)  the  algorithm  works  regardless 
of  whether  the  number  of  relay  nodes  in  the  network  is 
more  than  or  less  than  the  number  of  source-destination 
pairs;  (ii)  the  final  achievable  rate  for  each  source- 
destination  pair  is  guaranteed  to  be  no  less  than  that 
under  direct  transmissions;  (iii)  the  algorithm  is  able  to 
find  the  optimal  objective  regardless  of  initial  relay  node 
assignment. 
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•  We  provide  a  sketch  of  a  possible  implementation  of  the 
ORA  algorithm.  Some  practical  issues  and  overhead  in 
the  implementation  are  discussed. 

B.  Paper  Organization 

In  Section  II,  we  discuss  related  work  and  contrast  them 
with  this  paper.  Section  III  gives  a  brief  overview  of  coopera¬ 
tive  communications,  so  as  to  set  the  context  of  our  study.  In 
Section  IV,  we  describe  the  relay  node  assignment  problem 
in  a  cooperative  ad  hoc  network  environment.  Section  V 
presents  our  ORA  algorithm.  In  Section  VI,  we  give  a  proof 
of  optimality  for  ORA.  Section  VII  presents  numerical  results, 
and  Section  VIII  presents  a  sketch  of  how  ORA  can  be 
implemented.  Section  IX  concludes  this  paper. 

II.  Related  Work 

The  concept  of  cooperative  communications  can  be  traced 
back  to  the  three-terminal  communication  channel  (or  a  relay 
channel)  in  [20]  by  Van  Der  Meulen.  Shortly  after.  Cover  and 
El  Gamal  studied  the  general  relay  channel  and  established  an 
achievable  lower  bound  for  data  transmission  [4].  These  two 
seminal  works  laid  down  the  foundation  for  the  present-day 
research  on  cooperative  communications  that  can  be  broadly 
classified  into  the  following  three  categories. 

(a)  Physical  Layer  Schemes.  Current  research  on  CC  aims 
to  exploit  distributed  antennas  on  other  nodes  in  the  network. 
This  has  resulted  in  several  protocols  at  the  physical  layer  [5], 
[7],  [8],  [10],  [14],  [16],  [17].  These  protocols  describe  various 
ways  through  which  nodes  can  cooperate  at  the  physical  layer. 

In  [10],  Laneman  et  al.  studied  the  mutual  information 
between  a  pair  of  nodes  using  a  third  cooperating  node 
under  the  so-called  fixed  relaying  schemes  (AF  or  DF).  The 
underlying  physical  layer  model  for  CC  in  this  paper  is  based 
on  these  two  schemes.  In  addition  to  fixed  relaying  schemes, 
the  authors  also  presented  selection  relaying,  in  which  nodes 
can  switch  between  AF  or  DF  (depending  on  instantaneous 
channel  conditions),  and  incremental  relaying,  which  utilizes 
limited  feedback  from  the  receiving  node  to  further  improve 
the  performance  of  CC. 

In  [7],  Gunduz  and  Erkip  studied  an  opportunistic  coopera¬ 
tion  scheme  in  which  a  feedback  channel  among  cooperating 
nodes  can  be  used  to  share  channel  state  information  and 
help  perform  power  control.  The  authors  showed  that  the 
performance  of  DF  improves  when  power  control  is  employed. 
This  kind  of  opportunistic  DF  is  an  alternative  to  the  physical 
layer  fixed  DF  considered  in  our  work. 

Another  alternative  physical  layer  scheme  could  be  the 
delay-tolerant  DF  presented  in  [5].  In  this  delay-tolerant  DF 
scheme,  distributed  space-time  codes  are  used  to  address  the 
issue  of  asynchrony  (transmission  delay)  among  cooperative 
transmitters. 

Additionally,  in  [8]  and  [14],  authors  studied  multi-hop 
cooperative  protocols  that  involve  cooperation  among  multiple 
transmitting  nodes  along  the  path.  In  [16]  and  [17],  the 
authors  performed  an  in-depth  study  on  the  practical  issues 
of  implementing  user  cooperation  in  a  conventional  CDMA 
system. 


(b)  Network  Layer  Schemes  for  Multi-hop  Networks. 

Recent  efforts  on  CC  at  the  network  layer  include  [9],  [15], 
[22],  In  [9],  Khandani  et  al.  studied  minimum  energy  routing 
problem  (for  a  single  message)  by  exploiting  both  wireless 
broadcast  advantage  and  CC.  However,  their  proposed  solu¬ 
tions  cannot  provide  any  performance  guarantee  for  general 
ad  hoc  networks.  In  [22],  Yeh  and  Berry  aimed  to  generalize 
the  well  known  maximum  differential  backlog  policy  [18]  in 
the  context  of  CC.  They  formulated  a  challenging  nonlinear 
program  with  exponential  number  of  variables  that  character¬ 
izes  the  network  stability  region,  but  only  provided  solutions 
for  a  few  simple  network  topologies.  In  [15],  Scaglione  et  al. 
proposed  two  architectures  for  multi-hop  cooperative  wireless 
networks.  Under  these  architectures,  nodes  in  the  network 
can  form  multiple  cooperative  clusters.  They  showed  that 
the  network  connectivity  can  be  improved  by  using  such 
cooperative  clusters.  However,  problems  related  with  optimal 
routing  and  relay  node  assignment  were  not  discussed  in  their 
work. 

(c)  Relay  Node  Assignment  for  Ad  hoc  Networks.  The 

most  relevant  research  to  our  work  (i.e.  relay  node  assignment) 
include  [1],  [2],  [13],  [21],  [24],  In  [24],  Zhao  et  al.  showed 
that  for  a  single  source-destination  pair,  in  the  presence  of 
multiple  relay  nodes,  it  is  sufficient  to  choose  one  “best” 
relay  node,  instead  of  multiple  relay  nodes.  This  result  is 
interesting,  as  it  paves  the  way  for  research  on  assigning  no 
more  than  one  relay  node  to  a  source-destination  pair,  which 
is  the  setting  that  we  have  adopted  in  this  paper.  In  [21], 
Wang  et  al.  showed  how  game  theory  can  be  used  by  a 
single  session  to  select  the  best  cooperative  relay  node.  In  [1], 
Bletsas  et  al.  proposed  a  distributed  scheme  for  relay  node 
selection  based  on  the  instantaneous  channel  conditions  at  the 
relay  node.  In  contrast  to  [1],  [21],  and  [24],  our  paper  is 
not  limited  to  a  single-session,  and  considers  the  relay  node 
assignment  for  multiple  competing  sessions  with  the  goal  of 
maximizing  the  minimum  data  rate  among  all  of  them.  In  [13], 
Ng  and  Yu  studied  an  important  utility  maximization  problem 
for  the  joint  optimization  of  relay  node  selection,  cooperative 
communications,  and  resource  allocation  in  a  cellular  network. 
However,  their  solution  procedure  has  non-polynomial  running 
time.  In  [2],  Cai  et  al.  studied  relay  node  selection  and  power 
allocation  for  AF-based  wireless  relay  networks,  and  proposed 
a  heuristic  solution.  Additionally,  both  [2]  and  [13]  have 
different  objectives  from  our  work. 

III.  Cooperative  Communications:  A  Primer 
The  essence  of  cooperative  communications  is  best  ex¬ 
plained  by  a  three-node  example  in  Fig.  1.  In  this  figure,  node 
s  is  the  source  node,  node  d  is  the  destination  node,  and  node 
r  is  a  relay  node.  Transmission  from  s  to  d  is  done  on  a  frame- 
by-frame  basis.  Within  a  frame,  there  are  two  time  slots.  In 
the  first  time  slot,  source  node  s  makes  a  transmission  to  the 
destination  node  d.  Due  to  the  broadcast  nature  of  wireless 
communications,  this  transmission  is  also  overheard  by  the 
relay  node  r.  In  the  second  time  slot,  node  r  forwards  the 
data  received  in  the  first  time  slot  to  node  d.  Note  that  such  a 
two-slot  structure  is  necessary  for  cooperative  communications 
due  to  the  half-duplex  nature  of  most  wireless  transceivers. 
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Fig.  1.  A  three-node  schematic  for  cooperative  communication. 


In  this  section,  we  give  expressions  for  achievable  data  rate 
under  cooperative  communications  and  direct  transmissions 
(i.e.,  no  cooperation).  For  cooperative  communications,  we 
consider  both  amplify-and-forward  (AF)  and  decoded-and- 
forward  (DF)  modes  [10]. 

Amplify-and-Forward  (AF)  Under  this  mode,  let  hsd,  hsr, 
hrd  capture  the  effects  of  path-loss,  shadowing,  and  fading 
between  nodes  s  and  d,  s  and  r,  and  r  and  d,  respectively. 
Denote  Zd[  1]  and  Zd\ 2]  the  zero-mean  background  noise  at 
node  d  in  the  first  time  slot  and  second  time  slot,  respectively, 
both  with  variance  a2.  Denote  zr  [1]  the  zero-mean  background 
noise  at  node  r  in  the  first  time  slot,  with  variance  a2. 

Denote  xs  the  signal  transmitted  by  source  node  s  in  the 
first  time  slot.  Then  the  received  signal  at  destination  node  d, 
ysd ,  can  be  expressed  as 

Usd,  —  hsdXs  ~h  Z^[l]  ,  (1) 

and  the  received  signal  at  the  relay  node  r,  ysr,  is 

t/sr  =  hsrXs  T  ^r[l]  •  (2) 

In  the  second  time  slot,  relay  node  r  transmits  to  destination 
node  d.  The  received  signal  at  d,  yrd ,  can  be  expressed  as 

drd  —  ford  *  &r  *  Vsr  T  ^rf[2]  , 

where  ar  is  the  amplifying  factor  at  relay  node  r  and  ysr  is 
given  in  (2).  Thus,  we  have 

yrd  —  hrdOtr  '  ( hsrxs  -f-  zr[l])  T  £(z[2]  .  (3) 

The  amplifying  factor  ar  at  relay  node  r  should  satisfy  power 
constraint  a2(\hsr\2Pa  +  a2)  =  Pr,  where  Ps  and  Pr  are  the 
transmission  powers  at  nodes  s  and  r,  respectively.  So,  ar  is 
given  by 

2  _  Pr 

(Xr  ~  \hsr\2Ps  +  <j2  ‘ 

We  can  re-write  (1),  (2)  and  (3)  into  the  following  compact 
matrix  form 


Y  =  Ha  s  +  BZ  , 


where 


Y  = 

Usd 

,  H  = 

hsd 

Upd 

OLrhrdhSr 

0  10’ 

"  Mi] 

B  = 

(y.f'hf'ci  0  1 

,  and  Z  = 

*d[  1] 

.  zd[ 2]  . 

(4) 


It  has  been  shown  in  [10]  that  the  above  channel,  which 
combines  both  direct  path  (s  to  d)  and  relay  path  (s  to  r  to  d), 
can  be  modeled  as  a  one-input,  two-output  complex  Gaussian 


noise  channel.  The  achievable  data  rate  Caf(s,  r,  d)  from  s  to 
d  can  be  given  by 


CAF(s,r,d)  =  y  log2[det(I  +  (PsHHt)(BF;[ZZt]Bt)-1)]  , 

(5) 

where  W  is  the  bandwidth,  det(-)  is  the  determinant  function, 
I  is  the  identity  matrix,  the  superscript  “f”  represents  the 
complex  conjugate  transposition,  and  E[-]  is  the  expectation 
function. 

After  putting  (4)  into  (5)  and  performing  algebraic  ma¬ 


nipulations,  we  have  CAF{s,r,d)  =  —  log  (  1  +  ^-j\hsd\ 


Ps 


Ps\hsr\2  PrlKd? 


Ps<Xd\hsr\2 
\hsd\ 
have 


V 

D  2|,  |2  , — 2  2  i  •  Denote  SNRsd  = 

Prv2\hrd\2  +  a^ajj 

SNRsr  =  | i\hsr\ 2,  and  SNRrd  =  ^\hrd\2 ■  We 


Caf(s,  r,d)  =W  ■  IA P(SNRsd,  SNRsr,  SNRrd)  ,  (6) 

where  IAF( SNRsd,  SNRsr,  SNRrd)  =  \  log2  ^1  +  SNRsd 

,  SNR,,.SNRT.rf  ~) 

+  SNRsr+SNRrd+l )  ' 

Decode-and-Forward  (DF)  Under  this  mode,  relay  node  r 
decodes  and  estimates  the  received  signal  from  source  node  s 
in  the  first  time  slot,  and  then  transmits  the  estimated  data  to 
destination  node  d  in  the  second  time  slot.  The  achievable  data 
rate  for  DF  under  the  two  time-slot  structure  is  given  by  [10] 
as 


CDF(s,  r,d)=W-  ID  P(SNRsd)  SNRsr,  SNRrd )  ,  (7) 

where 

/DF(SNRsd,  SNRsr,  SNRrc()  =  ^  min{log2(l  +  SNRsr), 

log2(l  +  SNRsd  + SNRrd)}.  (8) 

Note  that  /ap(-)  and  /dp(-)  are  increasing  functions  of  Ps 
and  Pr,  respectively.  This  suggests  that,  in  order  to  achieve  the 
maximum  data  rate  under  either  mode,  both  source  node  and 
relay  node  should  transmit  at  maximum  power.  In  this  paper, 
we  let  Ps  =  Pr  =  P. 

Direct  Transmission  When  cooperative  communications 
(i.e.,  relay  node)  is  not  used,  source  node  s  transmits  to 
destination  node  d  in  both  time  slots.  The  achievable  data  rate 
from  node  s  to  node  d  is 


Cr,(s,d)  =  W  log2 (1  +  SNRsd)  . 

Based  on  the  above  results,  we  have  two  observations. 
First,  comparing  Cap  (or  Cop)  to  Cd,  it  is  hard  to  say  that 
cooperative  communications  is  always  better  than  the  direct 
transmission.  In  fact,  a  poor  choice  of  relay  node  could  make 
the  achievable  data  rate  under  cooperative  communications  to 
be  lower  than  that  under  direct  transmission.  This  fact  under¬ 
lines  the  significance  of  relay  node  selection  in  cooperative 
communications.  Second,  although  AF  and  DF  are  different 
mechanisms,  the  capacities  for  both  of  them  have  the  same 
form,  i.e.,  a  function  of  SNRsd,  SNRsr,  and  SNR,Y/.  Therefore, 
a  relay  node  assignment  algorithm  designed  for  AF  is  also 
applicable  for  DF.  In  this  paper,  we  develop  a  relay  node 
assignment  algorithm  for  both  AF  and  DF.  Table  I  lists  the 
notation  used  in  this  paper. 


(ft)  Potential  Relay  Node 
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TABLE  I 
Notation 


Symbol 

Definition 

Or  (s»,rj) 

Achievable  rate  for  Si~d{  pair  when  relay  node  rj  is 
used 

Or  (sj ,  0) 

Achievable  rate  for  Si~di  pair  under  direct  transmission 

r1  . 
wmin 

The  minimum  rate  among  all  source-destination 
pairs 

huv 

Effect  of  path-loss,  shadowing,  and  fading  from  node 
u  to  node  v 

Afs 

Set  of  source  nodes  in  the  network 

Afr 

Set  of  relay  nodes  in  the  network 

Ns 

=  | Afs  | ,  number  of  source  nodes  in  the  network 

Nr 

=  |  J\fr  | ,  number  of  relay  nodes  in  the  network 

N 

Number  of  all  the  nodes  in  the  network 

p 

Maximum  transmission  power 

r3 

The  j- th  relay  node,  rj  E  A fr 

R-ii>  (^i ) 

The  relay  node  assigned  to  Si  under  'ip 

Si 

The  i-th  source  node.  Si  E  Jfs 

O'!/;  {r  j  ) 

The  source  node  that  uses  rj  under  ip 

SNRU„ 

The  signal  noise  ratio  between  nodes  u  and  v 

W 

Channel  bandwidth 

Xs 

Signal  transmitted  by  node  s 

Vuv 

Received  signal  at  node  v  (form  node  u ) 

Zv[t ] 

Background  noise  at  node  v  during  time  slot  t 

Qj' 

Amplifying  factor  at  relay  r 

<*l 

Variance  of  background  noise  at  node  v 

A  solution  for  relay  node  assignment 

IV.  The  Relay  Node  Assignment  Problem 

Based  on  the  background  in  the  last  section,  we  consider 
relay  node  assignment  problem  in  a  network  setting.  There 
are  N  nodes  in  an  ad  hoc  network,  with  each  node  being 
either  a  source  node,  a  destination  node,  or  a  potential  relay 
node  (see  Fig.  2).  In  order  to  avoid  interference,  we  assume 
that  orthogonal  channels  are  available  in  the  network  (e.g., 
using  OFDMA),  which  is  proposed  for  cooperative  communi¬ 
cations  [10],  The  channel  gain  from  node  u  to  v  is  captured  by 
variable  huv.  Denote  Ms  =  (si,  S2,  •  •  • ,  sjv3  }  the  set  of  source 
nodes,  M,j  =  (di,  cfe,  ■  ■  • ,  d/vd}  the  set  of  destination  nodes, 
and  Mr  =  {ti,  r2,  •  •  • ,  rNr}  the  set  of  relays  (see  Fig.  2).  We 
consider  unicast  transmission  where  every  source  node  Sj  is 
paired  with  a  destination  node  di,  i.e.,  Nrj  =  Ns.  We  also 
consider  that  each  node  is  equipped  with  a  single  transceiver 
and  can  transmit/receive  within  one  channel  at  a  time.  We 
assume  that  each  node  can  only  serve  a  unique  role  of  source, 
destination,  or  relay.  That  is,  Nr  N  2 Ns.  Further,  we 
assume  that  a  session  utilizes  one  relay  node  for  CC  [24]. 

Note  that  a  source  node  may  not  always  get  a  relay  node. 
There  are  two  possible  scenarios  in  which  this  may  happen. 
First,  there  may  not  be  sufficient  number  of  relay  nodes  in  the 
network  (e.g.,  Nr  <  Ns).  In  this  case,  some  source  nodes  will 
not  have  relay  nodes.  Second,  even  if  there  are  enough  relay 
nodes,  a  sender  may  choose  not  to  use  a  relay  node  if  it  leads 
to  a  lower  data  rate  than  direct  transmission  (see  discussion  at 
the  end  of  Section  III). 

We  now  discuss  the  objective  function  of  our  problem. 
Although  different  objectives  can  be  used,  a  widely-used 
objective  for  CC  is  to  increase  the  achievable  data  rate  of 
individual  sessions.  For  the  multi-session  network  environment 
considered  in  this  paper  (see  Fig.  2),  each  source-destination 
pair  will  have  a  different  achievable  data  rate  after  we  apply 
a  relay  node  assignment  algorithm.  So,  a  plausible  objective 
is  to  maximize  the  minimum  data  rate  among  all  the  source- 


^  Sender  Q  Receiver 


Fig.  2.  A  cooperative  ad  hoc  network  consisting  of  source  nodes,  destination 
nodes,  and  relay  nodes. 

destination  pairs. 

More  formally,  denote  IZ(si)  the  relay  node  assigned  to  Sj, 
and  S(r:l  )  as  the  source  node  that  uses  r3 .  For  both  AF  and 
DF,  the  achievable  data  rate  of  the  session  can  be  written  as 
(see  Section  III) 

W/R(SNRSi)d4,  SNRSi!TC(s.),  SNRK(s.);d.)  , 

with  /r(-)  =  F af ( * )  when  AF  is  employed,  and  /r(-)  =  Idf(-) 
when  DF  is  employed.  In  case  Si  does  not  use  a  relay,  we 
denote  IZ(si)  =  0,  and  the  data  rate  is  the  achievable  rate 
under  direct  transmission,  i.e., 

Cr(s,,0)  =  CD(si,di )  ■ 

Combining  both  these  cases,  we  have 

cR{si,n(si))  = 

W/R(SNRSi)di,SNRSiiTC(si), 

SNRn{siUi)  if  7e(Si)^0  (9) 

Wlog(l  +  SNRSiidi)  ifn(si)  =  Q 

Note  that  we  do  not  list  di  in  function  Cr(s,;,  IZ(si))  since  for 
each  source  node  s,,  the  corresponding  destination  node  d,  is 
deterministic. 

Denote  Cmin  as  our  objective  function,  which  is  the  mini¬ 
mum  rate  among  all  source  nodes.  That  is, 

Cmin  —  min{CR(Si,  .  Si  (E  TVs}. 

Our  objective  is  to  find  an  optimal  relay  node  assignment  for 
all  the  source-destination  pairs  such  that  Cmin  is  maximized. 

In  subsequent  sections,  we  present  a  polynomial  time  so¬ 
lution  to  the  relay  node  assignment  problem  along  with  a 
correctness  proof. 

V.  An  Optimal  Relay  Assignment  Algorithm 
A.  Basic  Idea 

The  optimal  polynomial-time  algorithm  we  will  present  is 
called  Optimal  Relay  Assignment  (ORA)  algorithm.  Figure  3 
shows  the  flow  chart  of  the  ORA  algorithm. 

Initially,  the  ORA  algorithm  starts  with  a  random  but 
feasible  relay  node  assignment.  By  feasible,  we  mean  that 
each  source-destination  pair  can  be  assigned  at  most  one  relay 
node  and  that  a  relay  node  can  be  assigned  only  once.  Such 
initial  feasible  assignment  is  easy  to  construct,  e.g.,  direct 
transmission  between  each  source-destination  pair  (without  the 
use  of  a  relay)  is  a  special  case  of  feasible  assignment. 

Starting  with  this  initial  assignment,  ORA  adjusts  the  as¬ 
signment  during  each  iteration,  with  the  goal  of  increasing 
the  objective  function  Cmin.  Specifically,  during  each  iteration, 
ORA  identifies  the  source  node  that  corresponds  to  Cmin. 
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Fig.  4.  An  example  tree  topology  in  ORA  algorithm  for  finding  a  better 
solution. 


Fig.  3.  A  flow  chart  of  the  ORA  algorithm. 

Then,  ORA  helps  this  source  node  to  search  a  better  relay 
such  that  this  “bottleneck”  data  rate  can  be  increased.  In  the 
case  that  the  selected  relay  is  already  assigned  to  another 
source  node,  further  adjustment  of  relay  node  for  that  source 
node  is  necessary  (so  that  its  current  relay  can  be  released). 
Such  adjustment  may  have  a  chain  effect  on  a  number  of 
source  nodes  in  the  network.  It  is  important  that  for  any 
adjustment  made  on  a  relay  node,  the  affected  source  node 
should  still  maintain  a  data  rate  larger  than  Cmm .  There  are 
only  two  outcomes  from  such  search  in  an  iteration:  (i)  a  better 
assignment  is  found,  in  which  case,  ORA  moves  on  to  the  next 
iteration;  or  (ii)  a  better  assignment  cannot  be  found,  in  which 
case,  ORA  terminates. 

There  are  two  key  technical  challenges  we  aim  to  address 
in  the  design.  First,  for  any  non-optimal  solution,  the  algo¬ 
rithm  should  be  able  to  find  a  better  solution.  As  a  result, 
upon  termination,  the  final  assignment  is  optimal.  Second, 
its  running  time  must  be  polynomial.  We  will  show  that 


ORA  addresses  both  problems  successfully.  Specifically,  we 
show  the  complexity  of  the  ORA  algorithm  is  polynomial 
in  Section  V-D.  We  will  also  give  a  correctness  proof  of  its 
optimality  in  Section  VI. 

B.  Algorithm  Details 

In  the  beginning,  ORA  algorithm  performs  a  “preprocess¬ 
ing”  step.  In  this  step,  for  each  source-destination  pair,  the 
source  node  ,s,;  considers  each  relay  node  r:]  in  the  network 
and  computes  the  corresponding  data  rate  C^{si,rj)  by  (9). 
Each  source  node  ,s,  also  computes  the  rate  CR(Si,0)  by  (9) 
under  direct  transmissions  (i.e.,  without  the  use  of  a  relay 
node).  After  these  computations,  each  source  node  s*  can 
identify  those  relay  nodes  that  can  offer  an  increase  in  its 
data  rate  compared  to  direct  transmissions,  i.e.,  those  relays 
with  CR(Si,  Tj)  >  Cr(Si,  0).  Obviously,  it  only  makes  sense  to 
consider  these  relays  for  CC.  In  the  case  that  no  relay  can  offer 
any  increase  of  data  rate  compared  to  direct  transmissions,  we 
will  just  employ  direct  transmissions  for  these  source  nodes. 

After  the  preprocessing  step,  we  enter  the  initial  assignment 
step.  The  objective  of  this  step  is  to  obtain  an  initial  feasible 
solution  for  ORA  algorithm  so  that  it  can  start  its  iteration. 
In  the  pre-processing  step,  we  have  already  identified  the  list 
of  relay  nodes  for  each  source  node  that  can  increase  its  data 
rate  compared  to  direct  transmission.  We  can  randomly  assign 
a  relay  node  from  this  list  to  a  source  node.  Note  that  once  a 
relay  node  is  assigned  to  a  source  node,  it  cannot  be  assigned 
again  to  another  source  node.  Thus,  if  there  is  no  relay  node 
available  to  a  source  node,  then  this  source  node  will  simply 
employ  direct  transmission  as  its  initial  assignment.  Upon  the 
completion  of  this  assignment,  each  source  node  will  have  the 
data  rate  no  less  than  that  under  direct  transmission. 

The  next  step  in  the  ORA  algorithm  is  to  find  a  better 
assignment,  which  represents  an  iteration  process.  This  is  the 
key  step  in  the  ORA  algorithm.  The  detail  of  this  step  is  shown 
in  the  bottom  portion  of  Fig.  3.  As  a  starting  point  of  this  step, 
ORA  algorithm  identifies  the  smallest  data  rate  C'rnm  among 
all  sources.  ORA  algorithm  aims  to  increase  this  minimum 
rate  for  the  corresponding  source  node,  while  having  all  other 
source  nodes  maintain  their  data  rates  above  Cm Without 
loss  of  generality,  we  use  Fig.  4  to  illustrate  a  search  process. 

•  Suppose  ORA  identifies  that  si  has  the  smallest  rate  Cjnm 
under  the  current  assignment  (with  relay  node  ri).  Then 
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Si  examines  other  relays  with  a  rate  larger  than  G'nlm.  If  it 
cannot  find  such  a  relay,  then  no  better  solution  is  found 
and  the  ORA  algorithm  terminates. 

In  case  of  a  tie,  i.e.,  when  two  or  more  source  nodes  have 
the  same  smallest  data  rate,  the  tie  is  broken  by  choosing 
the  source  node  with  the  highest  node  index. 

•  Otherwise,  i.e.,  if  there  are  better  relays,  we  consider 
these  relays  in  the  non-increasing  order  in  terms  of  data 
rate  (should  it  be  assigned  to  si).  That  is,  we  try  the  relay 
that  can  offer  the  maximum  possible  increase  in  data  rate 
first.  In  case  of  a  tie,  i.e.,  when  two  or  more  relay  nodes 
offer  the  same  maximum  data  rate,  the  tie  is  broken  by 
choosing  the  relay  node  with  the  highest  node  index. 

•  Suppose  that  source  node  Si  considers  relay  node  r2.  If 
this  relay  node  is  not  yet  assigned  to  any  other  source 
node,  then  r2  can  be  immediately  assigned  to  si.  In  this 
simple  case,  we  find  a  better  solution  and  the  current 
iteration  is  completed. 

•  Otherwise,  i.e.,  r2  is  already  assigned  to  a  source  node, 
say  52,  we  mark  r2  to  indicate  that  r2  is  “under  consid¬ 
eration"  and  check  whether  r2  can  be  released  by  s2. 

•  To  release  r2,  source  node  s2  needs  to  find  another  relay 
(or  use  direct  transmission)  while  making  sure  that  such 
new  assignment  still  has  its  data  rate  larger  than  C'mm . 
This  process  is  identical  to  what  we  have  done  for  si, 
with  the  only  (but  important)  difference  that  s2  will  not 
consider  a  relay  that  has  already  been  “marked”,  as  that 
relay  node  has  already  been  considered  by  a  source  node 
encountered  earlier  in  the  search  process  of  this  iteration. 

•  Suppose  that  source  node  s2  now  considers  relay  ?’3.  If 
this  relay  node  is  not  yet  assigned  to  any  source  node, 
then  r3  can  be  assigned  to  s2;  r2  can  be  assigned  to  si; 
and  the  current  iteration  is  completed.  Moreover,  if  the 
relay  under  consideration  by  s2  is  the  one  that  is  being 
used  by  the  source  node  that  initiated  the  iteration,  i.e., 
relay  r\,  then  it  is  easy  to  see  that  n  can  be  taken  away 
from  si.  A  better  solution,  where  r\  is  assigned  to  s2, 
and  r2  is  assigned  to  si,  is  found  and  the  current  iteration 
is  completed.  Otherwise,  we  mark  r 3  and  check  further 
to  see  whether  r3  can  be  released  by  its  corresponding 
source  node,  say  S3.  We  also  note  that  s2  can  consider 
direct  transmission  if  it  offers  a  data  rate  larger  than  C'mm . 

•  Suppose  that  S3  cannot  find  any  “unmarked"  relay  that 
offers  a  data  rate  larger  than  C'nlm,  and  its  data  rate  under 
direct  transmission  is  no  more  than  Cmin-  Then  s2  cannot 
use  ?’3  as  its  relay. 

•  If  any  “unmarked”  relay  that  offers  a  data  rate  larger 
than  Cmin  cannot  be  assigned  to  s2,  then  si  cannot  use 
r2  and  will  move  on  to  consider  the  next  relay  on  its 
non-increasing  rate  list,  say  r±. 

•  The  search  continues,  with  relay  nodes  being  marked 
along  the  way,  until  a  better  solution  is  found  or  no  better 
solution  can  be  found.  For  example,  in  Fig.  4,  Sq  finds 
a  new  relay  ry .  As  a  result,  we  have  a  new  assignment, 
where  rr  is  assigned  to  sq\  tq  is  assigned  to  S4;  and  r 4 
is  assigned  to  si. 

Note  that  the  “mark"  on  a  relay  node  will  not  be  cleared 


Main  algorithm 

1.  Perform  preprocessing  and  an  initial  relay  node  assignment. 

2.  Set  all  the  relay  nodes  in  the  network  as  "unmarked”. 

3.  Denote  s b  the  source  node  with  Cmin,*  the  smallest  data 
rate  among  all  source  nodes.  The  corresponding  destination 
node  of  Sf,  is  db  and  the  corresponding  relay  node  is  TZ(sb). 

4.  Find_Another_Relay  (S5,  TZ(sb),  Cm;n). 

5.  If  si,  finds  a  better  relay,  then  go  to  line  2. 

6.  Otherwise,  the  algorithm  terminates. 

Subroutines 

Find_Another_Relay(<S  (rj  ),rj,  Cj^): 

7.  For  every  "unmarked”  relay  r/.  with  Oj^(S(rj),  r^)  >  C'lnln . 
do  the  following  in  the  non-increasing  order  of  Cj^(S(rj),  rj.).* 

8.  Run  Check_Relay_Availability(ffc ,  Cm;n). 

9.  If  Tf.  is  available,  then  do  the  following: 

10.  Remove  relay  node  rj’s  assignment  to  <S(ry); 

11.  Assign  relay  node  r/,.  to  S(rj). 

12.  Otherwise,  continue  on  to  next  rj.  and  go  to  line  8. 

13.  If  all  relays  are  unavailable,  then  <S ( Vj )  cannot  find  another  relay. 

Check_Relay_Availability(rj ,  Cmjn): 

14.  If  r-j  is  not  assigned  to  any  source  node,  then  Vj  is  available. 

15.  If  r-j  =  TZ(sb)  or  rj  =  0.  then  r-f  is  available. 

16.  Otherwise, 

17.  Set  r-j  as  “marked”. 

18.  Run  Find_Another_Relay  (S(vj),Tj,Cmjn). 

19.  If  S(rj)  can  find  another  relay,  then  r-j  is  available. 

20.  Otherwise  r-j  is  unavailable. 

*  A  tie  is  broken  by  choosing  the  node  with  the  largest  node  index. 

Fig.  5.  Pseudocode  for  the  ORA  algorithm. 

throughout  the  search  process  in  the  same  iteration.  We  call 
this  the  “linear  marking”  mechanism.  These  marks  will  only 
be  cleared  when  the  current  iteration  terminates  and  before 
the  start  of  the  next  iteration.  A  pseudocode  for  the  ORA 
algorithm  is  shown  in  Fig.  5. 

We  now  use  an  example  to  illustrate  the  operation  of  the 
ORA  algorithm,  in  particular,  its  “linear  marking”  mechanism. 
Readers  who  already  understood  the  ORA  algorithm  can  skip 
this  example. 

Example  1:  Suppose  that  there  are  seven  source-destination 
pairs  and  seven  relay  nodes  in  the  network. 

Table  11(a)  shows  the  data  rate  for  each  source  node  Sj 
when  relay  node  r-j  is  assigned  to  it.  The  symbol  0  indicates 
direct  transmission.  Also  shown  in  Table  11(a)  is  an  initial 
relay  node  assignment,  which  is  indicated  by  an  underscore 
on  the  intersecting  row  (sj)  and  column  (rj).  Note  that  the 
preprocessing  step  before  the  initial  assignment  ensures  that 
the  data  rate  for  each  source-destination  pair  in  the  initial 
assignment  is  no  less  than  that  under  direct  transmission. 

Under  the  initial  relay  node  assignment  in  Table  11(a),  source 

53  is  identified  as  the  bottleneck  source  node  Sb  with  the 
smallest  rate  of  C'Tnm  =  13.  Since  consideration  of  relay  nodes 
is  performed  in  the  order  of  non-increasing  (from  largest  to 
smallest)  data  rate  for  the  source  node  under  consideration, 
r4  is  therefore  considered  for  S3.  Since  r4  is  already  assigned 
to  source  node  s2,  we  “mark”  r4  now.  Now  s2  needs  to  find 
another  relay.  But  any  other  relay  (or  direct  transmission)  will 
result  in  a  data  rate  no  greater  than  the  current  objective  value 
Cniin  =  13.  This  means  that  r4  cannot  be  taken  away  from  s2. 
Since  r 4  does  not  work  out  for  S3,  S3  will  then  consider  the 
next  relay  node  that  offers  the  second  largest  data  rate  value, 
i.e.,  relay  node  r 7.  Since  r 7  is  already  assigned  to  sender  S4, 
we  “mark”  r-  now.  Next,  ORA  algorithm  will  check  to  see  if 

54  can  find  another  relay.  It  turns  out  that  none  of  the  relay 
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TABLE  II 

An  example. 


(a)  Initial  relay  node  assignment. 


0 

r  1 

T2 

r3 

r4 

r5 

r7 

Si 

14 

7 

24 

5 

14 

15 

17 

9 

S2 

9 

8 

10 

11 

20 

10 

12 

11 

— >  S3 

11 

10 

13 

17 

21 

8 

9 

19 

S4 

12 

8 

9 

12 

11 

10 

9 

18 

S5 

10 

9 

18 

19 

24 

9 

13 

23 

S6 

7 

18 

12 

6 

11 

11 

17 

20 

S7 

16 

1 

9 

4 

14 

19 

8 

12 

(b)  Assignment  after  the  first  iteration. 


0 

ri 

T2 

V3 

7*4 

T5 

VQ 

7*7 

->•  si 

14 

7 

24 

5 
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(c)  Assignment  after  the  second  iteration. 
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(d)  Final  assignment  upon  termination. 
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nodes  except  r7  can  offer  a  data  rate  larger  than  the  current 
Cmin  to  s 4.  As  a  result,  r^  cannot  be  taken  away  from  S4. 
Source  node  S3  will  now  check  for  the  relay  node  that  offers 
next  largest  rate,  i.e.,  r 3.  Since  r 3  is  already  assigned  to  sender 
S5,  we  “mark”  r 3  now.  Next,  ORA  algorithm  checks  to  see  if 
S5  can  find  another  relay.  Then  S5  checks  relay  nodes  in  non¬ 
increasing  order  of  data  rate  values.  Since  r 4  (with  largest 
rate),  r7  (with  the  second  largest  rate),  and  r$  (with  the  third 
largest  rate)  are  all  marked,  they  will  not  be  considered.  The 
relay  with  the  fourth  largest  rate  is  r2,  which  offers  a  rate  of 
18  >  Cm in  =  13.  Moreover,  ?’2  is  the  relay  node  assigned  to 
Sb  =  S3-  Thus,  S5  can  choose  r-> ■  The  new  assignment  after 
the  first  iteration  is  shown  in  Table  11(b).  Now  the  objective 
value,  Cmin,  is  updated  to  15,  which  corresponds  to  si.  Before 
the  second  iteration,  all  markings  done  in  the  first  iteration  are 
cleared. 

In  the  second  iteration,  ORA  algorithm  will  identify  si  as 
the  source  node  with  a  minimum  data  rate  in  the  network. 
The  algorithm  will  then  perform  a  new  search  for  a  better 
relay  node  for  source  si.  Similar  to  the  first  iteration,  the 
assignments  for  other  source  nodes  can  change  during  this 
search  process,  but  all  assignments  should  result  in  data  rates 
larger  than  15. 

The  iteration  continues  and  the  final  assignment  upon  ter¬ 
mination  of  ORA  algorithm  is  shown  in  Table  11(d),  with  the 


optimal  (maximum)  value  of  Cmin  being  17.  ■ 

It  should  be  clear  that  ORA  works  regardless  of  whether 
Nr  >  Ns  or  Nr  <  Ns.  For  the  latter  case,  i.e.,  the  number  of 
relay  nodes  in  the  network  is  less  than  the  number  of  source 
nodes,  it  is  only  necessary  to  consider  relay  node  assignment 
for  a  reduced  subset  of  Nr  source  nodes,  where  the  data  rate 
of  each  source  in  this  subset  under  direct  transmission  is  less 
than  the  data  rate  of  those  (Ns  —  Nr)  source  nodes  not  in  this 
subset.  As  a  result,  in  the  case  of  Ns  >  Nr,  ORA  will  run 
even  faster  due  to  a  smaller  problem  size. 

C.  Caveat  on  the  Proposed  Marking  Mechanism 

We  now  re-visit  the  marking  mechanism  in  the  ORA 
algorithm.  Although  different  marking  mechanisms  may  be 
designed  to  achieve  the  optimal  objective,  the  algorithm 
complexity  under  different  marking  mechanisms  may  differ 
significantly.  In  this  section,  we  first  present  a  marking  mech¬ 
anism,  which  appears  to  be  a  natural  approach  but  leads  to  an 
exponential  complexity  for  each  iteration.  Then  we  discuss  our 
proposed  marking  mechanism  and  show  its  linear  complexity 
for  each  iteration. 

A  natural  approach  is  to  perform  both  marking  and  un¬ 
marking  within  an  iteration.  This  approach  is  best  explained 
with  an  example.  Again,  let’s  look  at  Fig.  4.  Source  node 
si  first  considers  r-2 .  Since  r2  is  being  considered  by  si  in 
the  new  solution  and  is  used  by  S2  in  the  current  solution, 
7"2  is  marked.  Source  node  S2  considers  r%,  which  is  already 
assigned  to  S3.  Since  S3  cannot  release  r 3  without  reducing 
its  data  rate  below  the  current  Cmin,  this  branch  of  search  is 
futile  and  si  now  considers  a  different  relay  node  r,\.  Since 
r4  is  currently  assigned  to  S4,  we  mark  r,\  and  try  to  find  a 
new  relay  for  S4.  Now  the  question  is:  shall  we  remove  those 
marks  on  r  >  and  r 3  that  we  put  on  earlier  in  the  process  within 
this  iteration?  Under  this  natural  approach,  r2  and  r 3  should 
be  unmarked  so  that  they  can  be  considered  as  candidate 
relay  nodes  for  S4  in  its  search.  Similarly,  when  we  try  to 
find  a  relay  for  sq,  relay  nodes  r-2 ,  t-j  ,  r4  and  r-,  should  be 
unmarked  so  that  they  can  be  considered  as  candidate  relay 
nodes  for  Sq,  in  addition  to  r7.  It  is  not  hard  to  show  that 
such  marking/unmarking  mechanism  will  consider  all  possible 
assignments  and  can  guarantee  to  find  an  optimal  solution 
upon  termination.  However,  the  complexity  of  such  approach 
is  exponential  within  each  iteration. 

In  contrast,  under  the  ORA  algorithm,  there  is  no  unmarking 
mechanism  within  an  iteration.  That  is,  relay  nodes  that  are 
marked  earlier  in  the  search  process  by  some  source  nodes  will 
remain  marked.  As  a  result,  any  relay  node  will  be  considered 
at  most  once  in  the  search  process,  which  leads  to  a  linear 
complexity  for  each  iteration.  Unmarking  for  all  nodes  is 
performed  only  at  the  end  of  an  iteration  so  that  there  is  a 
clean  start  for  the  next  iteration. 

An  immediate  question  regarding  our  marking  mechanism 
is:  how  could  such  a  “linear  marking”  lead  to  an  optimal 
solution,  as  it  appears  that  many  possible  assignments  that  may 
increase  Cmin  are  not  considered.  This  is  precisely  the  question 
that  we  will  address  in  Section  VI,  where  we  will  prove  that 
ORA  can  guarantee  that  its  final  solution  is  optimal. 


D.  Complexity  Analysis 

We  now  analyze  the  computational  complexity  of  ORA 
algorithm.  During  each  iteration,  due  to  the  “linear  marking” 
mechanism  in  our  algorithm,  a  relay  node  is  checked  for 
its  availability  at  most  once.  Thus,  the  complexity  of  each 
iteration  is  0{Nr). 

Now  we  examine  the  maximum  number  of  iterations  that 
ORA  can  execute.  The  number  of  improvements  in  data  rate 
that  an  individual  source  node  can  have  is  limited  by  Nr. 
As  a  result,  in  worst  case,  the  number  of  iterations  that 
the  algorithm  can  go  through  are  0(NsNr).  This  makes  the 
overall  complexity  of  ORA  algorithm  to  be  0(NsN 0). 

VI.  Proof  of  Optimality 

In  this  section,  we  give  a  correctness  proof  of  the  ORA 
algorithm.  That  is,  upon  the  termination  of  the  ORA  algorithm, 
the  solution  (i.e.,  objective  value  and  the  corresponding  relay 
node  assignment)  is  optimal. 

Our  proof  is  based  on  contradiction.  Denote  0  the  final 
solution  obtained  by  the  ORA  algorithm,  with  the  objective 
value  being  G'nnn.  For  0,  denote  the  relay  node  assigned  to 
source  node  s,  as  7 Z^(si).  Conversely,  for  0,  denote  the  source 
node  that  uses  relay  node  r,  as  S^(rj). 

We  now  assume  that  there  exists  a  solution  0  better  than  0. 
That  is,  the  objective  value  by  0,  denoted  as  C'Tnm,  is  greater 
than  that  by  0,  i.e.,  Cmin  >  Cm in.  For  -0,  we  denote  the  relay 
node  assigned  to  source  node  ,s,  as  Conversely,  for  0, 

we  denote  the  source  node  that  uses  relay  node  ry  as  Sy  (ry ) . 

The  key  idea  in  the  proof  is  to  exploit  the  marking  status 
of  relay  nodes  at  the  end  of  its  last  iteration,  which  is  a  non¬ 
improving  iteration.  Specifically,  in  the  beginning  of  this  last 
iteration,  ORA  will  select  a  “bottleneck”  source  node,  which 
we  denote  as  Sb-  ORA  will  then  try  to  improve  the  solution 
by  searching  for  a  better  relay  node  for  this  bottleneck  source 
node.  Since  the  last  iteration  is  a  non-improving  iteration, 
ORA  will  not  find  a  better  solution,  and  thus  will  terminate. 
We  will  show  that  IZy(sb)  is  not  marked  at  the  end  of  the  last 
iteration  of  ORA.  On  the  other  hand,  by  assuming  that  there 
exists  a  better  solution  0  than  0,  we  will  show  that  IZy(sb) 
will  be  marked  at  the  end  of  the  last  iteration  of  ORA.  This 
leads  to  a  contradiction  and  thus  0  cannot  exist.  We  begin  our 
proof  with  the  following  fact. 

Fact  1:  For  the  bottleneck  source  node  Sb  under  0,  its  relay 
node  lZy>(sb)  is  not  marked  at  the  end  of  the  last  iteration  of 
the  ORA  algorithm. 

Proof:  In  the  ORA  algorithm,  a  relay  node  r:i  is  marked 
only  if  r-j  0  IZy(sb)  (see  Check_Relay_Availability()  in 
Fig.  5).  Thus,  IZy(sb)  cannot  be  marked  at  the  end  of  the 
last  iteration  of  the  ORA  algorithm.  ■ 

Fact  1  will  be  the  basis  for  contradiction  in  our  proof  for 
Theorem  1,  the  main  result  of  this  section. 

Now  we  present  the  following  three  claims,  which  recur¬ 
sively  examine  relay  node  assignment  under  0.  First,  for  the 
relay  node  assigned  to  Sb  in  0,  i.e.,  7 Zy(sb),  we  have  the 
following  claim. 

Claim  1:  Relay  node  7 Zy(sb)  must  be  marked  at  the  end  of 
the  last  iteration  of  the  ORA  algorithm.  Further,  it  cannot  be 
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Fig.  6.  The  sequence  of  nodes  under  analysis  in  the  proof  of  optimality. 


0  and  must  be  assigned  to  some  source  node  under  solution 

0- 

Proof:  Since  0  is  a  better  solution  than  0,  we  have 
CR(sb,TZy(sb))  >  ©mi,,  >  C'min.  Thus,  by  construction,  ORA 
will  consider  the  relay  node  IZyfsbY s  availability  for  Sb  in 
its  last  iteration.  Since  ORA  algorithm  cannot  find  a  better 
solution  in  its  last  iteration,  relay  IZy(sb)  should  be  marked 
and  then  the  outcome  for  checking  TZyfsbfs  availability  must 
be  unavailable.  By  “linear  marking”,  the  mark  on  TZy(sb) 
will  not  be  cleared  throughout  the  search  process  in  the  last 
iteration.  Thus,  the  relay  node  TZy(sb)  is  marked  at  the  end 
of  the  last  iteration  of  ORA  algorithm. 

We  now  prove  the  second  statement  by  contradiction.  If 
TZy(sb)  is  0,  then  Sb  will  choose  0  in  the  last  iteration  since 
it  can  offer  CR(sb,7Zy(sb))  >  Cmin-  But  this  contradicts  to 
the  fact  that  we  are  now  in  the  last  iteration  of  ORA,  which 
is  a  non-improving  iteration.  So  TZy(sb)  cannot  be  0.  Further, 
since  we  proved  that  TZy(sb)  is  marked  at  the  end  of  the  last 
iteration  of  the  ORA  algorithm,  it  must  be  assigned  to  some 
source  node  already.  ■ 

By  the  definition  of  Sy{-),  we  have  that  TZy(sb)  is  assigned 
to  source  node  Sy(lZy(sb))  in  solution  0.  To  simplify  nota¬ 
tion,  define  function  Qy(-)  as 

M-)  =  <*v(  7M-))-  (10) 

Thus,  relay  node  IZy(sb)  is  assigned  to  source  node  Gy{sb) 
in  0  (see  top  portion  of  Fig.  6). 

Since  IZy(sb)  0  7 Zy(sb),  they  are  assigned  to  different 
source  nodes  in  0,  i.e.,  Gy>{sb)  0  Sb-  Now,  we  recursively 
investigate  the  relay  node  assigned  to  source  Gy(sb)  under 
solution  0,  i.e.,  7 Zy(Gy(sb)).  We  have  the  following  claim 
(also  see  Fig.  6). 

Claim  2:  Relay  node  TZy(Gy(sb))  must  be  marked  at  the 
end  of  the  last  iteration  of  the  ORA  algorithm.  Further,  it 
cannot  be  0  and  must  be  assigned  to  some  source  node  under 
solution  0. 

The  proofs  for  both  statements  in  this  claim  follow  the  same 
token  as  that  for  Claim  1. 
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Again,  by  the  definition  of  <S^(-),  we  have  that  relay  node 
is  assigned  to  source  node  S^(Uj(G^(sb)))  in 
solution  ip.  By  (10),  we  have  source  (TZ^ [Gy, (sb)))  = 

Gii>(Gii>(sb))-  To  simplify  the  notation,  we  define  function 

Gil)  as 

Gl(-)=G^(GA  •))• 

Thus,  relay  node  7?.^  (<?,/, (sf,))  is  assigned  to  source  node 
G2,(sb)  in  ip-  Now  we  have  two  cases:  source  node  G2,(sb)  may 
or  may  not  be  a  node  in  {sb,  G^(sb)}.  If  source  node  G2,(sb) 
is  a  node  in  {sb,  Gij>(sb)},  then  we  terminate  our  recursive 
procedure.  Otherwise,  we  can  further  consider  its  relay  node 
in  ip. 

In  general  we  can  use  the  following  notation. 

Gi/j{Sb)  =  Sb, 

Gl{sb)=  GAG^isb))  (k>  l).  (ll) 

Since  the  numbers  of  source  nodes  are  finite,  our  recursive 
procedure  will  terminate  in  finite  steps.  Suppose  that  we 
terminate  at  k  =  n. 

Following  the  same  token  for  Claims  1  and  2,  we  can 
obtain  a  similar  claim  for  each  of  the  relay  nodes  7£^(<?^,(sb))> 
Kj(Gl(sb)),  •  •  -,n4xg^sb)),  •  •  -,n^(sb))  (see  Fig.  6). 
Thus,  we  can  generalize  the  statements  in  Claims  1  and  2  for 
relay  node  lZ^(G^(sb))  and  have  the  following  claim. 

Claim  3:  Relay  node  T^-^{G^{sb))  must  be  marked  at  the 
end  of  the  last  iteration  of  the  ORA  algorithm.  Further,  it 
cannot  be  0  and  must  be  assigned  to  some  source  node  under 
solution  ip,  k  =  0, 1,  2,  •  •  • ,  n. 

Proof:  Since  'ip  is  a  better  solution  than  ip,  we  can  say  that 
Cr(<4(s&),7^(^(s6)))  ^  l^min  >  Cmin-  Note  that  G^{sb)  is 
some  source  node  in  the  solution  ip  obtained  by  ORA,  whereas 
7l,ij,(Glp1(sb))  is  the  relay  node  assigned  to  this  source  node  in 
the  hypothesized  better  solution  ip.  Our  goal  is  to  show  that 
ORA  should  have  marked  this  relay  node  in  its  last  iteration. 

Since  CR{G^{sb),lZ^{G^{sb)))  >  Cmin  and  Tlj(G§(sb))  is 
not  assigned  to  QP(sb)  in  the  last  iteration  of  ORA,  then  by 
construction  of  ORA,  ORA  must  have  checked  7^(f/^(s&))’s 
availability  for  Q%  (sb)  during  the  last  iteration,  then  marked  it, 
and  then  determined  it  to  be  unavailable  for  Q\{sb).  Moreover, 
due  to  “linear  marking”,  this  mark  on  IZ^iG^jfsb))  should  be 
there  after  the  last  iteration  of  ORA.  Thus,  we  can  conclude 
that  IZ^fQ^fsb))  is  marked  at  the  end  of  the  last  iteration  of 
the  ORA  algorithm. 

We  now  prove  the  second  statement  by  contradiction.  If 
TZ^iG^f  sb))  is  0,  then  Q^(sb)  will  choose  0  in  the  last  iteration 
since  it  can  offer  CR(G§(sb),  H^(G^,(sb))  >  Cmin,  and  finally 
sb  will  be  able  to  get  a  better  relay  node.  But  this  contradicts 
with  the  fact  that  this  last  iteration  is  a  non-improving  iteration. 
So,  'R.,ijj(G,j>{sb))  cannot  be  0.  Further,  since  we  proved  that 
'R'AG^(sb))  is  marked  at  the  end  of  the  last  iteration  of  the 
ORA  algorithm,  it  must  be  assigned  to  some  source  node 
already.  ■ 

Referring  to  Fig.  6,  we  have  Claim  3  for  a  set  of  relay  nodes 
TZ^{sb),n^{Gy,{sb)),  ■  ■  Tl^G^isb))-  Our  recursive  proce¬ 
dure  terminates  at  'Tl^{G7^{sb))  because  its  assigned  source 
node  in  solution  ip  is  a  node  in  {sb,  G,jAsb)  > ' ' ' ,  GPp{sb)}.  We 


are  now  ready  to  prove  the  following  theorem,  which  is  the 
main  result  of  this  section. 

Theorem  1:  Upon  the  termination  of  the  ORA  algorithm, 
the  obtained  solution  ip  is  optimal. 

Proof:  Under  Claim  3,  we  proved  that  the  relay  node 
IZ.ifGPf  sb))  is  assigned  to  some  source  node  in  solution  ip 
obtained  by  ORA.  Since  our  recursive  procedure  terminates 
at  7£^((/^;(sb))>  its  assigned  source  node  in  solution  ip  is 
a  node  in  {sb,Gy,{sb),-  ■  ■  ,Glp{sb)}.  But  we  also  know  that 
under  ip,  source  nodes  Gyj{sb),  G^,{sb),  G^{sb),  ■  • G^sf) 
have  relay  nodes  K^(sb),  ^(G^{sb)),  ^(0$(«b)).  ■  • 
TZ^G^isb)),  respectively.  Thus,  fc${G${8b))  is  the  only 
relay  node  that  can  be  assigned  to  sb  in  solution  ip.  On  the 
other  hand,  relay  node  assigned  to  sb  in  solution  ip  is  denoted 
by  IZ^(sb).  Thus,  we  have  H^(G^{sb))  =  7 Z^(sb). 

Now,  Claim  3  states  that  7 ^(G^(sb))  must  be  marked  after 
the  last  iteration,  whereas  Fact  1  states  that  the  relay  node 
assigned  to  the  bottleneck  source  node,  i.e.,  7 Z^(sb),  cannot 
be  marked.  Since  both  lZ^(sb)  and  7 ^(f?^(sb))  are  the  same 
relay  node,  we  have  a  contradiction.  Thus  our  assumption  that 
there  exists  a  solution  ip  better  than  ip  does  not  hold.  The  proof 
is  complete.  ■ 

Note  that  the  proof  of  Theorem  1  does  not  depend  on  the 
initial  assignment  in  ORA.  So  we  have  the  following  important 
property. 

Corollary  1.1:  Under  any  feasible  initial  relay  node  assign¬ 
ment,  the  ORA  algorithm  can  find  an  optimal  relay  node 
assignment. 

VII.  Numerical  Results 

In  this  section,  we  present  some  numerical  results  to  demon¬ 
strate  the  properties  of  the  ORA  algorithm. 

A.  Simulation  Setting 

We  consider  a  100-node  cooperative  ad  hoc  network.  The 
location  of  each  node  is  given  in  Table  III.  For  this  network,  we 
consider  both  the  cases  of  Nr  >  Ns  and  Nr  <  Ns.  In  the  first 
case,  we  have  30  source-destination  pairs  and  40  relay  nodes. 
While  in  the  second  case,  we  have  40  source-destination  pairs 
and  only  20  relay  nodes.  The  role  of  each  node  (either  as  a 
source,  destination,  or  relay)  for  each  case  is  shown  in  Figs.  7 
and  9,  respectively,  with  details  given  in  Table  III. 

For  the  simulations,  we  assume  W  =  10  MHz  bandwidth 
for  each  channel.  The  maximum  transmission  power  at  each 
node  is  set  to  1  W.  Each  relay  node  employs  AF  for  cooper¬ 
ative  communications.  We  assume  that  hsd  only  includes  the 
path  loss  component  between  nodes  s  and  d  and  is  given  by 
\hsd\2  =  ||s  —  d||~4,  where  ||s  —  d||  is  the  distance  (in  meters) 
between  these  two  nodes  and  4  is  the  path  loss  index.  Note 
that  the  working  of  the  ORA  algorithm  does  not  depend  on  the 
mode  of  CC  and  the  channel  gain  model.  As  long  as  channel 
gains  and  achievable  rates  are  known,  ORA  will  give  optimal 
assignment.  For  the  AWGN  channel,  we  assume  the  variance 
of  noise  is  10~10  W  at  all  nodes. 
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TABLE  III 

Locations  and  roles  of  all  the  nodes  in  the  network. 


Node 

Location 

Role 

Node 

Location 

Role 

Node 

Location 

Role 

Case  1 

Case  2 

Case  1 

Case  2 

Case  1 

Case  2 

(75,  500) 

si 

si 

(220,  190) 

d<± 

<^4 

(380,  370) 

r*7 

S31 

(170,430) 

S2 

S2 

(660,  190) 

d5 

d$ 

(300,  350) 

78 

7  8 

(170, 500) 

S3 

S3 

(430,  630) 

de 

dft 

(410,650) 

rg 

S33 

(250,  650) 

S4 

S4 

(180,620) 

d'j 

dy 

(470,  500) 

rio 

d40 

(400, 550) 

(750,  625) 

d8 

ds 

(660,  525) 

rn 

«39 

(340, 230) 

S6 

S6 

(310,  480) 

dg 

dg 

(600,  425) 

ri2 

S40 

(390, 150) 

S7 

S7 

(1100, 180) 

dio 

di  o 

(510,  200) 

ri3 

S38 

(460, 280) 

S8 

S8 

(1110,360) 

dn 

dn 

(575, 325) 

ri4 

7*14 

(700,  500) 

S9 

S9 

(875,  600) 

di2 

di2 

(750,  560) 

ri5 

7*15 

(750,  360) 

Sio 

Sio 

(700,  300) 

di3 

di3 

(800,  360) 

riQ 

7-16 

(800,  90) 

s  11 

sn 

(650,  550) 

di4 

di4 

(860,  260) 

ri7 

7*17 

(900, 160) 

S 12 

S12 

(740,  170) 

d  15 

dis 

(980,  450) 

ris 

7*18 

(1125,300) 

$13 

S13 

(410,810) 

di6 

^16 

(950,310) 

r\9 

7-19 

(1000,  340) 

S 14 

314 

(550,  1100) 

di7 

di7 

(950,  200) 

r20 

d37 

(1025, 540) 

S15 

S15 

(150,  790) 

dis 

dis 

(100, 1000) 

T21 

s  32 

(100, 1120) 

S16 

316 

(210,  1110) 

di9 

dig 

(310,980) 

T22 

7-12 

(150,920) 

S17 

S17 

(530,  720) 

dg  0 

dgg 

(250,  800) 

T23 

d32 

(330, 1110) 

S18 

S18 

(800,  1140) 

^21 

d2l 

(460, 1010) 

T24 

7*13 

(450, 890) 

S19 

319 

(1080, 1100) 

d22 

dg2 

(610,930) 

r25 

^34 

(650, 1050) 

S20 

S20 

(940,  790) 

d23 

d23 

(680,  760) 

T26 

S34 

(700, 640) 

S21 

S21 

(1360,640) 

d24 

d24 

(700,  900) 

T27 

7*20 

(820,  880) 

$22 

322 

(1280, 1120) 

d25 

d25 

(910, 1120) 

?-28 

d35 

(1150, 1060) 

S23 

S23 

(1260,350) 

d26 

d26 

(970,  970) 

7-29 

S35 

(1480, 1120) 

S24 

324 

(1500,50) 

d27 

d27 

(1360,910) 

7  30 

79 

(1160,  720) 

S25 

325 

(1450,605) 

d28 

d28 

(1200,920) 

7-31 

7*11 

(1050,  50) 

S26 

S26 

(1030,910) 

dgg 

dgg 

(1250,690) 

7-32 

d36 

(1350,  450) 

S 27 

327 

(1150,230) 

d30 

d30 

(1290, 180) 

733 

7*10 

(1380,  110) 

S28 

S28 

(80, 370) 

ri 

d3l 

(150,360) 

7-34 

7-5 

(1500,  800) 

S29 

S29 

(110,280) 

7*2 

T2 

(1380,380) 

7-35 

7-7 

(1500,  300) 

S30 

330 

(160,300) 

73 

r  3 

(1220,60) 

736 

s37 

(200,  50) 

di 

di 

(280,  520) 

7*4 

V4 

(1190,510) 

7-37 

S36 

(520, 240) 

d2 

d2 

(375,  580) 

7*5 

d39 

(500,  40) 

7-38 

d38 

(40,  100) 

d3 

d3 

(385,  450) 

76 

(50,805) 

739 

d33 

(1510,920) 

7-40 

7*1 

0  Senders  Q  Receivers  (Q)  Potential  Relays 


Fig.  7.  Topology  for  a  100-node  network  for  Case  1  (Nr  >  Ns),  with 
Ns  =  30  and  Nr  =  40. 


B.  Results 

Case  1:  Nr  >  Ns.  In  this  case  (see  Fig.  7),  we  have  30 
source-destination  pairs  and  40  relay  nodes. 

Under  ORA,  after  preprocessing,  we  start  with  an  initial 
relay  node  assignment  in  the  first  iteration.  Such  initial  as¬ 
signment  is  not  unique.  But  regardless  of  the  initial  relay 
node  assignment,  we  expect  the  objective  value  to  converge 
to  the  optimum  (by  Corollary  1.1).  To  validate  this  result,  in 
Table  IV,  we  show  the  results  of  running  the  ORA  algorithm 
under  two  different  initial  relay  node  assignments,  denoted  as 


Iteration  Number 

Fig.  8.  Case  1  (Nr  >  Ns):  The  objective  value  Cmjn  at  each  iteration  of 
ORA  algorithm  under  two  different  initial  relay  node  assignments. 

I  and  II  (see  Table  IV). 

In  Table  IV,  the  second  column  shows  the  data  rate  for 
each  source-destination  pair  under  direct  transmissions.  Note 
that  the  minimum  rate  among  all  pairs  is  1.83  Mbps,  which 
is  associated  with  S7.  The  third  to  fifth  columns  are  results 
under  initial  relay  node  assignment  I  and  sixth  to  eighth 
columns  are  results  under  initial  relay  node  assignment  II. 
The  symbol  0  denotes  direct  transmissions.  Note  that  initial 
relay  node  assignments  I  and  II  are  different.  As  a  result,  the 
final  assignment  is  different  under  I  and  II.  However,  the  final 
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TABLE  IV 

Optimal  assignments  for  Case  1  (Nr  >  Na)  under  two  different 

INITIAL  RELAY  NODE  ASSIGNMENTS. 


•  Senders  Q  Receive,s  ©  Potential  Relays 


Fig.  9.  Topology  for  a  100-node  network  for  Case  2  (Nr  <  Ns),  with 
Ns  =  40  and  Nr  =  20. 


objective  value  (i.e.,  C'mm )  under  I  and  II  is  identical  (4.43 
Mbps). 

Figure  8  shows  the  objective  value  Cjrml  at  each  iteration 
under  initial  relay  node  assignments  I  and  II.  Under  either 
initial  relay  node  assignments  I  or  II,  C'mm  is  a  non-decreasing 
function  of  iteration  number.  Note  that  a  higher  initial  value 
of  Cm in  does  not  mean  that  ORA  will  converge  faster.  The 
increase  of  Crnm  by  cooperative  communications  over  direct 
transmissions  is  significant  (from  1.83  Mbps  to  4.43  Mbps). 


TABLE  V 

Optimal  assignments  for  Case  2  ( Nr  <  N3 )  under  two  different 

INITIAL  RELAY  NODE  ASSIGNMENTS. 


Fig.  10.  Case  2  (Nr  <  Ns):  The  objective  value  C'mjn  at  each  iteration  of 
ORA  algorithm  under  two  different  initial  node  assignments. 


source-destination  pairs  and  20  relay  nodes. 

Table  V  shows  the  results  of  this  case  under  two  different 


Case  2:  Nr  <  Ns. 


In  this  case  (see  Fig.  9),  we  have  40 


initial  relay  node  assignments  I  and  II.  The  second  column 
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TABLE  VI 

An  example  illustrating  the  importance  of  preprocessing. 


Sender 

Cd 

(Mbps) 

Without  Preprocessing  | 

Initial 

Final 

Final 

Rate 

(Mbps) 

si 

2.62 

rg 

7*3 

6.54 

S2 

4.60 

0 

0 

4.60 

S3 

3.81 

0 

r2 

8.73 

S4 

2.75 

rg 

7*4 

4.66 

S5 

3.15 

r  14 

7*14 

6.47 

S6 

4.17 

0 

7*6 

9.25 

S7 

1.83 

T6 

7*8 

4.76 

S8 

2.99 

0 

7*12 

7.22 

S9 

4.92 

0 

0 

4.92 

sio 

4.80 

0 

0 

4.80 

Sll 

4.13 

0 

7*20 

9.13 

S12 

3.24 

0 

ri8 

5.55 

S13 

3.68 

0 

7*17 

7.32 

S14 

4.23 

0 

ri6 

7.87 

S15 

2.62 

0 

r  19 

5.84 

Sl6 

3.30 

0 

7*22 

7.30 

S17 

4.17 

0 

7*24 

5.62 

S18 

6.03 

T23 

V21 

7.37 

Sl9 

8.76 

T39 

r39 

4.81 

S20 

6.95 

T26 

7*26 

7.25 

S21 

1.90 

0 

7*27 

4.90 

S22 

7.65 

r28 

7*28 

8.71 

S23 

7.55 

T29 

7*29 

11.26 

S24 

2.12 

0 

7*40 

4.43 

S25 

3.91 

0 

7*30 

5.87 

S26 

6.08 

T33 

7*33 

7.55 

S27 

3.61 

0 

7*34 

5.45 

S28 

2.04 

0 

7*35 

5.29 

S29 

2.33 

0 

7*31 

4.68 

S30 

6.60 

0 

0 

6.60 

in  Table  V  lists  the  data  rate  under  direct  transmissions.  As 
discussed  at  the  end  of  Section  V-B,  for  the  case  of  Nr  < 
Ns,  it  is  only  necessary  to  consider  relay  node  assignment 
for  Nr  =  20  source  nodes  corresponding  to  the  20  smallest 
achievable  rates  under  direct  transmission. 

Again  in  Table  V,  the  objective  value  C'm;n  is  identical  (3.80 
Mbps)  regardless  of  different  initial  relay  node  assignments 
(I  and  II).  Note  that  despite  the  difference  in  final  relay 
node  assignments  under  I  and  II,  the  objective  value  Cmm  is 
identical.  The  increase  of  C'Tnm  by  cooperative  communications 
over  direct  transmissions  is  significant  (from  1.83  Mbps  to 
3.80  Mbps). 

Figure  10  shows  the  objective  value  Cm ;n  at  each  iteration 
under  initial  relay  node  assignments  I  and  II.  Again,  we 
observe  that  in  Fig.  10,  Cm ;n  is  a  non-decreasing  function  of 
iteration  number  under  both  initial  relay  node  assignments  I 
and  II. 

Significance  of  Preprocessing:  Now  we  use  a  set  of 

numerical  results  to  show  the  significance  of  preprocessing 
in  our  ORA  algorithm.  We  consider  the  same  network  in 
Fig.  7  with  30  source-destination  pairs  and  40  relay  nodes. 
Now  we  remove  the  preprocessing  step  in  the  ORA  algorithm. 
As  an  example,  the  third  column  of  Table  VI  shows  an  initial 
assignment  without  first  going  through  the  preprocessing  step. 
Although  the  objective  value  Cm jn  also  reaches  the  same 
optimal  value  (4.43  Mbps)  as  that  in  Table  IV,  the  final  data 
rate  for  some  non-bottleneck  source  nodes  could  be  worse 
than  direct  transmissions.  For  example,  for  sig,  its  final  rate  is 


4.81  Mbps,  which  is  less  than  its  direct  transmission  rate  (8.76 
Mbps).  Such  event  is  undetectable  without  the  preprocessing 
step,  as  4.81  Mbps  is  still  greater  than  the  optimal  objective 
value  (4.43  Mbps). 

On  the  other  hand,  when  the  preprocessing  step  is  employed, 
ORA  can  ensure  that  the  final  rate  for  each  source-destination 
pair  is  no  less  than  that  under  direct  transmission,  as  shown 
in  Table  IV. 

VIII.  A  Sketch  of  a  Possible  Implementation 

In  this  section,  we  present  a  sketch  of  a  possible  implemen¬ 
tation  of  the  ORA  algorithm.  This  implementation  follows  the 
link-state  approach.  Note  that  although  a  link-state  approach 
is  not  considered  fully  distributed,  it  is  nevertheless  a  viable 
implementation,  as  evidenced  by  the  widespread  deployment 
of  OSPF  [12]  in  the  Internet  and  acceptance  of  OLSR  [3]  in 
wireless  ad  hoc  networks. 

A.  Ensuring  Identical  Optimal  Solution  at  Source  Nodes 

In  the  presentation  of  the  ORA  algorithm  in  Section  V,  we 
have  learned  that  the  ORA  algorithm  can  start  with  any  random 
initial  relay  node  assignment  and  can  still  obtain  an  optimal 
solution.  However,  in  the  link-state  based  implementation, 
each  source  node  in  the  network  will  run  ORA  independently 
on  its  own.  As  such,  the  randomness  in  initial  relay  assignment 
must  be  removed  in  implementation  so  as  to  ensure  that 
each  source  node  can  obtain  an  identical  optimal  solution. 
Otherwise,  we  may  run  into  a  situation  that  the  same  relay 
node  may  be  assigned  to  multiple  source  nodes. 

A  simple  way  to  ensure  identical  initial  relay  node  assign¬ 
ment  is  to  have  each  source  node  choose  direct  transmission, 
i.e.  0  as  its  initial  relay  assignment.  Given  such  identical  initial 
assignment  and  that  ORA  is  a  deterministic  algorithm,  each 
source  node  will  obtain  an  identical  final  optimal  solution. 

B.  Some  Implementation  Details 

Under  such  implementation,  each  relay  node  collects  its 
link  state  information  with  its  neighboring  source  nodes;  each 
destination  node  also  collects  its  link  state  information  with 
its  source  node  and  neighboring  relay  nodes.  To  do  this, 
each  source  node  sends  a  broadcast  packet  to  its  neighboring 
relay  nodes  and  its  destination  node;  each  relay  node  sends 
a  broadcast  packet  to  its  neighboring  destination  nodes.  As 
shown  in  [6]  by  Gollakota  and  Katabi,  this  broadcast  packet 
transmission  can  be  used  by  the  receiver  of  each  wireless 
link  to  accurately  determine  the  link  state.  It  was  also  shown 
in  [6]  that  such  an  approach  is  practically  feasible,  and  was 
demonstrated  in  their  implementation  of  802.11  receivers. 
However,  we  point  out  that  in  an  uncontrolled  environment, 
estimating  channel  state  is  not  trivial. 

Upon  obtaining  the  link-state  information,  each  relay  and 
destination  node  will  distribute  such  information  to  all  the 
source  nodes  in  the  network.  This  will  ensure  that  each  source 
node  will  have  global  link-state  information.  Such  link-state 
dissemination  can  be  achieved  by  using  one  of  the  many  effi¬ 
cient  flooding  techniques  (see  e.g.  [23])  for  wireless  networks. 
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The  overhead  of  this  operation  is  small  when  compared  to 
potential  gain  in  achievable  rate  in  optimal  assignment  (see 
Section  VIII-C). 

Once  each  source  node  has  global  link-state  information,  it 
can  now  run  ORA  locally,  with  an  identical  initial  assignment 
as  discussed  in  Section  VIII-A.  As  discussed,  the  final  optimal 
solution  obtained  at  each  source  node  will  be  identical. 


C.  Overhead 


An  important  consideration  in  our  implementation  is  the 
overhead  incurred  in  distributing  link-state  information  in  the 
network.  This  can  be  measured  by  comparing  such  overhead 
with  the  potential  gain  in  achievable  rate  in  the  optimal 
solution.  We  now  analyze  such  overhead  and  show  that  the 
ratio  between  the  two  is  small,  thus  affirming  the  efficacy  of 
our  proposed  implementation. 

First  of  all,  Mandke  et  al.  [11]  conducted  extensive  exper¬ 
iments  and  showed  that  the  CSI  for  the  10  MHz  band  in  2.4 
GHz  spectrum  changes  every  300  mSec  on  average.  Their 
experiments  showed  that  CSI  distribution  does  not  need  to 
be  performed  more  frequently  than  every  300  mSec  (or  0.3 
Sec). 

To  run  the  ORA  algorithm  locally,  each  source  node  must 
obtain  global  link-state  information.  This  can  be  done  by 
having  all  relay  nodes  and  destination  nodes  flood  their  local 
link-state  information  in  the  network.  To  estimate  an  upper 
bound  for  such  flooding  overhead,  we  assume  that  32  bits 
(commonly  used  for  floating  variables)  are  used  to  represent 
each  link-state  value.  Then  the  total  link-state  information 
collected  at  each  relay  node  has  Or  =  Ns  x  32  bits.  Similarly, 
the  total  link-state  information  collected  at  each  destination 
node  has  Od  =  (Nr  +  1)  x  32  bits.  As  a  result,  the  total 
overhead  (in  b/s)  due  to  flooding  at  every  node  is: 


O  = 


Nr  ■  Or  +  Nd  ■  Od 
CL3 


(12) 


As  an  example,  for  Case  1  in  the  numerical  results  in  Sec¬ 
tion  VII-B,  it  can  be  shown  that  the  total  overhead  is  170.83 
Kb/s  at  each  node.  On  the  other  hand,  the  gain  in  the  bottleneck 
data  rate  by  ORA  is  4.43  —  1.83  =  2.6  Mb/s.  The  ratio  between 
the  two  is  only  6.4%.  That  is,  the  overhead  is  much  less  than 
the  gain  of  CC. 

We  acknowledge  that  in  some  environments,  the  overhead 
could  be  large  if  CSI  in  the  network  varies  on  a  smaller  time 
scale.  Under  such  environment,  fast  and  efficient  dissemination 
of  CSI  remains  an  open  problem. 


IX.  Conclusion 

Cooperative  communications  is  a  powerful  communication 
paradigm  to  achieve  spatial  diversity.  However,  the  perfor¬ 
mance  of  such  communication  paradigm  hinges  upon  the 
assignment  of  relay  nodes  in  the  network.  In  this  paper, 
we  studied  this  problem  in  a  cooperative  ad  hoc  network 
environment,  where  multiple  source-destination  pairs  compete 
for  the  same  pool  of  relay  nodes.  Our  objective  is  to  assign  the 
available  relay  nodes  to  different  source-destination  pairs  so  as 
to  maximize  the  minimum  data  rate  among  all  the  pairs.  The 


main  contribution  of  this  paper  is  a  polynomial  time  optimal 
algorithm  that  achieves  this  objective.  A  novel  idea  in  this 
algorithm  is  a  “linear  marking”  mechanism,  which  is  able  to 
achieve  linear  complexity  at  each  iteration.  We  gave  a  formal 
proof  of  optimality  for  the  algorithm  and  used  numerical 
results  to  demonstrate  its  efficacy. 

Although  we  offered  a  sketch  of  a  possible  implementation 
of  ORA,  a  number  of  issues  remain  challenging  in  practice. 
In  particular,  fast  and  efficient  method  for  collecting  and 
disseminating  CSI  in  moderate  and  large  sized  networks 
remain  an  open  problem.  Nevertheless,  the  theoretical  results 
presented  here  can  be  used  as  a  performance  benchmark  for 
other  proposed  solutions  in  practice. 
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